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Within the framework of quantization of the macroscopic electromagnetic field, a master equation 
describing both the resonant dipole-dipole interaction (RDDI) and the resonant atom-field interac- 
tion (RAFI) in the presence of dispersing and absorbing macroscopic bodies is derived, with the 
relevant couplings being expressed in terms of the surroundings-assisted Green tensor. It is shown 
that under certain conditions the RDDI can be regarded as being governed by an effective Hamil- 
tonian. The theory, which applies to both weak and strong atom-field coupling, is used to study 
the resonant energy exchange between two (two-level) atoms sharing initially a single excitation. 
In particular, it is shown that in the regime of weak atom-field coupling there is a time window, 
where the energy transfer follows a transfer-rate law of the type obtained by ordinary second-order 
perturbation theory. Finally, the spectrum of the light emitted during the energy transfer is studied 
and the line splittings are discussed. 

PACS numbers: 42.50.Ct, 42.50.Fx, 42.60.Da, 82.20.Rp 



I. INTRODUCTION 



It is well known that an initially excited molecule can 
irreversibly transfer its energy to unexcited species lo- 
cated nearby - a fundamental process involved in ultra- 
fast photochemistry, excitation transport in Langmuir- 
Blodgett films, determination of molecular structures, 
etc.. For (acceptor) molecules that exhibit a sufficiently 
broad (quasi-)continuum of vibronic states, the energy 
transfer may be regarded as being irreversible, and the 
transfer rate can be calculated by using standard second- 
order perturbation theory (see, e.g., [Q and references 
therein). If instead of molecules, the donor and acceptor 
species are (two-level) atoms, the transfer can be reversed 
in principle, leading to a back-and-forth energy exchange 
between the atoms [^|. This process can be thoroughly 
described by a nonperturbative treatment. When the in- 
teratomic interaction is weak, however, the probability 
of the return transfer can be so small that the energy 
exchange is essentially again a one-way process. In ad- 
dition to that, the atoms always feature some natural 
linewidths, which can be regarded as continua of states. 
An aim of the present article is to establish the connec- 
tion between the intermolecular energy transfer and the 
spontaneous-decay-assisted interatomic energy exchange. 

The underlying mechanism of intermolecular and in- 
teratomic energy transfer can be both RDDI and RAFI 
H . In free space and for transfer distances R below the 
wavelength of real-photon emission, the RDDI described 
in terms of an instantaneous Coulomb interaction poten- 
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tial with the familiar R~ 3 distance law is typically con- 
sidered as the leading mechanism. For transfer distances 
exceeding the wavelength of real-photon emission, the 
RDDI may be thought of as resulting from an extended 
off-resonant atom-field interaction, and the effect of the 
RAFI needs to be considered carefully. In particular, 
in the multipolar-coupling quantum electrodynamics [Q , 
the near-zone and far-zone results simply arise as limits 
of a unified theory based on the full atom-field coupling 
(in electric-dipole approximation). 

Similar to the single-atom spontaneous decay, the in- 
teratomic energy exchange can be strongly affected by 
the presence of material surroundings. While RDDI and 
RAFI in free space have been more or less well under- 
stood, there have been a number of open problems when 
material surroundings are present. In particular, it is 
not clear whether the near-zone R~ 3 law is still generally 
valid. Another problem that needs investigation is the 
regime of strong RAFI, i.e., when suitably chosen mate- 
rial surroundings form a high-Q (micro-)resonator, which 
gives rise to a strong atom-field coupling accompanied by 
quantum Rabi oscillations. These and related problems 
have recently been very actual in view of proposals to 



entangle atoms (see, e.g.,_ 
turn logic gates (see, e.g. 



and to implement quan- 
7| ) , using atoms in a cavity, 



impurities in the condensed phase, or quantum dots in a 
semiconductor . 

For two atoms that interact with cavity-type modes, 
a Hamiltonian for treating both RDDI and RAFI was 
suggested in . In the Hamiltonian, an effective dipole- 
dipole interaction energy, which may be thought of as re- 
sulting from a perturbative treatment of the off-resonant 
atom-field coupling, is added to the exact on-resonant 
atom-field interaction energy. Unfortunately, the rel- 
evant coupling strengths are introduced as parameters 
rather than physically determined system quantities. 
Moreover, the underlying mode expansion concept makes 
the theory inapplicable for the realistic case of absorbing 



2 



and dispersing material surroundings. Inclusion of mate- 
rial absorption and dispersion is especially important in 
designing quantum computing systems where dissipation 
and dccoherence have to be strictly kept under control. 

In what follows we give a more rigorous approach to 
the problem, basing on a recently developed quantiza- 
tion scheme for the electromagnetic field in the presence 
of (linear) absorbing and dispersing macroscopic bodies 
(see |j) and references therein). In this scheme, an op- 
erator noise polarization source is introduced into the 
Maxwell equations, and the electromagnetic field oper- 
ators are given by a source-quantity representation in 
terms of the classical Green tensor and a continuous 
set of bosonic vector fields, which play the role of the 
fundamental variables of the composite system of the 
electromagnetic field plus the media. The Green ten- 
sor, which can be calculated in a purely classical man- 
ner, is essentially determined by the (Kramers-Kronig 
consistent) space- and frequency-dependent complex per- 
mittivity, which characterizes the relevant properties of 
the macroscopic bodies. The medium-assisted electro- 
magnetic field can then be coupled to additional atomic 
systems, and the minimal-coupling Hamiltonian can be 
constructed, from which the corresponding multipolar- 
coupling Hamiltonian can be derived by means of a uni- 
tary transformation. 

By using the multipolar-coupling Hamiltonian in 
electric-dipole approximation, in Section II an equation 
of motion for the reduced density operator of a system 
that consists of two-level atoms and that part of the 
medium-assisted electromagnetic field that resonantly in- 
teracts with the atoms is derived. It is shown that 
when the (shifted) atomic transition frequencies differ by 
amounts small compared to the scale of the frequency 
variation of the Green tensor, the system can be regarded 
as being governed by a Hamiltonian, with the RDDI be- 
ing described by a well-defined effective interaction en- 
ergy. The problem of resonant energy exchange between 
two two- level atoms is studied in Section III, with special 
emphasis on the rate regime. The spectrum of the light 
emitted during the energy exchange is considered in Sec- 
tion [fv|. Finally, some concluding remarks are given in 
Section [v|. 
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Here, f (r, u>) and V (r, ui) are bosonic fields which play the 
role of the fundamental variables of the electromagnetic 
field and the medium, including a reservoir necessarily as- 
sociated with the losses in the medium, G(r, r', uj) is the 
classical Green tensor, and e(r, uj) — £r(i", uj) + iei(r, uj) is 
the complex permittivity. There are no direct Coulomb 
forces between the particles in the multipolar-coupling 
Hamiltonian (|l|); all interactions are mediated by the 
medium-assisted electromagnetic field. 

With regard to the interaction of the atoms with the 
electromagnetic field, it is convenient to decompose the 
latter into an on-resonant part (denoted by J °° duo . . .) 

and an off-resonant part (denoted by/ °°dijj . . .). Let us 
consider the temporal evolution of the system that con- 
sists of the atoms and the on-resonant part of the elec- 
tromagnetic field. For this purpose we formally solve the 
Heisenberg equation of motion for f(r,o;) [and F(r,w)] 
and insert the result into the equation of motion for an 
arbitrary system operator. Applying a (coarse-grained) 
Markov approximation to the slowly varying atomic vari- 
ables in the off-resonant frequency integrals, and assum- 
ing that the off-resonant free-field is initially (t — 0) pre- 
pared in the vacuum state, after some algebra, we derive 
the following equation of motion for the system density 
operator in the Schrodinger picture: 
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II. REDUCED DENSITY OPERATOR 
EQUATION 



A system of N two-level atoms [positions r^, tran- 
sition frequencies uja, transition dipole moments 
(A = 1, 2, N)] interacting with the electromagnetic 
field via electric-dipole transitions in the presence of dis- 
persing and absorbing bodies can be described by the 
multipolar-coupling Hamiltonian [Q, ^| 

H= [d 3 r[ dwfiwf f (r,w)f(r,w) + Vifc A ^ 
J Jo , 



where 
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= /d 3 r/ duhuP(r,u})f(r,u)+^±HujAVAz 
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dw 



d A E(r AjW )+H.c. , (5) 



and the notation Y^'a A' indicates that A'. In Eqs. 
(^) and (||) , the shifted atomic transition frequencies and 
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the resonant dipole-dipole coupling strengths are defined 
according to 



<^A = U A ~ <5 'A*4 



5a* a = S 
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(6) 
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S aa> -ZiT" / duj — " TT^r. (8) 



ttHsq J q c 2 uj — (+)Q A i 



[P - principal value]. The notation A* (A 1 *) means that 
d A (d^') in Eq. (||) has to be replaced with its complex 
conjugate (d^,). In Eq. (Q), both the RDDI and the 
RAFI are taken into account, without any restriction to 
the strength of the latter one. Applying the Kramers- 
Kronig relation to the Green tensor, we find from Eq. (||) 
that the relations 
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are (approximately) valid. 

From Eq. ([h]) it is seen that when the transition fre- 
quencies uiA and oja' of two atoms A and A' are different 
from each other, then |5>i»a'| and |#A'*a| are not sym- 
metric with respect to A and A' . Only if the differences 
l&A — Co A >\ are small compared with the frequency scale 
of variation of the Green tensor, so that 



J A*A> — V A'A 



St, 



(11) 



holds, this asymmetry can be disregarded. Then Eq. (Q) 
reduces to 



(jls - ^ ftbA*A>cr\vA> } 'J 
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Recalling the definition of Hs, we see that the motion of 
the system is governed by the effective Hamiltonian 

H cB = Jd 3 rJ dwRwft(r,w)f(r,w) 

+ ^ \Wb A a Az - ^ h5A*A>v\vA' 

A A,A> 

-J2 dw d A E(r j4 ,c J )+H.c. , (13) 



with E(r A ,u>) being given by Eq. (||). The Hamiltonian 
( |l3| ) is remarkable in several respects. Firstly, it applies 



to atoms surrounded by arbitrarily configured, disper- 
sive and absorptive media. Secondly, it goes beyond 
the rotating-wave approximation. Thirdly, it contains 
a RDDI energy that is explicitly expressed in terms of 
the medium-assisted Green tensor, according to Eq. ([h]). 
Note that both the single-atom transition frequency shift 
and the interatom RDDI result from the off-resonant 
atom-field coupling. 

Strong RDDI may be expected if the atoms are suffi- 
ciently close to each other. To give a simple example of 
the effect of material absorption, let us assume that the 
atoms are embedded in bulk material of complex permit- 
tivity s(lo). Using the bulk-material Green tensor (see, 
e.g., §), fromEq. ([[(]) we find in the short-distance limit 
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_ (14) 
(R,=r A ~r A >). In free space, Eq. (14) reduces to the well- 
known result that the RDDI simply corresponds to the 
(near-field) Coulomb-type interaction. Note that though 
the characteristic i? -3 distance dependence observed in 
free space is not changed by the bulk medium, it may be 
changed by surroundings with more complex geometry. 
In the long-distance limit, we find that 
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[n(ui) =nn(uj) + in\{uj) = ye(tJ)], i.e., the (harmonically 
modulated) R^ 1 dependence observed in free space 
is changed to an exponential decrease according to 
exp[— ni(Q A ')uj A iR/c] due to material absorption. Thus, 
material absorption can drastically reduce the RDDI 
strength with increasing mutual distance of the atoms. 



III. RESONANT ENERGY EXCHANGE BY A 
TWO-ATOM SYSTEM 

Let us assume that initially the atoms share a sin- 
gle excitation while the field is in the vacuum state. In 
the Schrodinger picture we may write, on omitting off- 
resonant terms, the state vector of the system in the form 
of 



■i(CuA —uj)t 



\UaW)) 



d 3 r / dujC Li (r,uj,t)e 



-i(u-w)t\T\ if 



|L)//(r,o;)|{0}) (16) 



(u> = \ J2a &a)- Here, \U A ) is the atomic state with the 
Ath atom in the upper state and all the other atoms in 
the lower state, and \L) is the atomic state with all atoms 
in the lower state. Accordingly, |{0}) is the vacuum state 
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of the rest of the system, and (r, u>)|{0}) is a state, 
where a single quantum is excited. 

From the effective Hamiltonian given in Eq. jl3|), we 
obtain the following system of coupled integrodifferential 
equations for the (slowly-varying) upper-state probabil- 
ity amplitudes C' A : 



When the two atoms are identical and have equiva- 
lent positions and dipole orientations with respect to the 
material surroundings such that the relations 

= r s . B , (25) 

Ta*b=T b *a, Sa'B=$B'A (26) 



C A (t) = iSA-A'e^-^CA'it) 

A' 
A'^A 

+ V] f dt 1 [°°duK A * A ,(t,t';u)C A ,{t'), (17) p a(B) (t) = \ [cosh (r 
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are valid, we have K-a* b —ICb*a, and Ta*b, T b *a, 5 a *b, 
and 5b* a are real quantities due to the reciprocity of the 
Green tensor. Then Eqs. ( |l9| ) and ( pp| ) yield 



,-r. 



where 



K AA ,(t,t';uj) 



1 



x d A lm G(r A , r A > , u)d A 



(18) 



We see that the RAFI is determined by the imaginary 
part of the Green tensor. In order to get insight into 
the atomic motion on the basis of closed solutions, let 
us consider two atoms and restrict our attention to the 
limiting cases of weak and strong RAFI. 



A. Weak atom-field coupling 

In the weak coupling regime, the integral expression 
in Eqs. ( |Tt| ) can be treated in a (coarse-grained) Markov 
approximation. We now make the simplifying assump- 
tion that the transition frequencies of the two atoms 
are nearly equal to each other, Cj a ~ Cjb , but allow for 
^A'AT^^B'B- The latter may happen, e.g., when the 
atoms A and B have different dipole matrix elements 
and/or different dipole orientations. It is then not diffi- 
cult to solve Eqs. (p"7j). In particular, if the atom B is 
initially in the lower state, Cs(i = 0) = 0, we obtain 

CA = ^{[-k{T A *A-T B * B ) + D] e D ^ 2 

+ [\{T A « A -T B , B ) + D\e D - t / 2 } 1 (19) 



(20) 



where 



^AA' 



d A lmG(r A ,r A >,oj A >)d A > , (21) 



K. AB = -\Y AB +i8 AB , (22) 
D = [\{T A , A - T B , B ) 2 + AK A * B K B * A ) 1/2 , (23) 



D± = -\(T A , A +T B , B )±D. 



(24) 



(t)=\C A( B)(t)\ 2 }. Eqs. (^ 



A'Bt) + (-)cos {25 A * B t)]e 

(24) [and thus 



Eq. ([271)1 are valid for arbitrary dispersing and absorb- 
ing material surroundings of the atoms. In particular, 
substituting in Eqs. ( ^l| ) and ( ^2|) for the Green tensor 
the vacuum Green tensor, Eq. (|27|) reduces to that one 
obtained in Ref. Q. Accordingly, using the Green ten- 
sor for absorbing bulk material, the result in Ref. |l(| is 
recognized. 

From Eq. (|27|) a damped oscillatory excitation ex- 
change between the two atoms is seen. For suf- 
ficiently small times, r B * B t<g^l, and strong RDDI, 
|<5a*s|^ F b * b , the oscillatory behavior dominates. This 
can typically be observed when the atoms are sufficiently 
near to each other, but can also be realized for more mod- 
erate distances with the interatom coupling being medi- 
ated by a high-Q medium-assisted field resonance [ [Tl| . 
In the opposite limit of weak resonant dipole-dipole cou- 
pling, P A (t) decreases monotonously while P B (t) features 
one peak which separates the regime of energy transfer 
from atom A to atom B at early times and the subse- 
quent decay of the excited state of atom B. 



B. Rate regime 

As already mentioned, for weak RDDI the energy 
transfer is one-way; that is, from atom A to atom B. 
In this transfer rate w\ can be defined according 

to 



W\ = 



dP B (t) 



dt 



where to is determined from the conditions that 
d 2 P B (t) _ n dP B (t) 



dH 



= 0, 



dt 



> 0. 



(28) 



(29) 



From Eq. (|2C|) [together with Eqs. (|2J) and (|2J) for 
negligibly small K A * B and K B - A therein] it then follows 
that 



Wi 



\Kb*a\ 2 D _ tn 
D 2 
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2Dt 



(D + + D^)e ut "\ ,(30) 
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1 —{(D + +D-f 
2D[{D+ + D-) 2 + AD+D-] 1 ' 2 }^. (31) 



Let us analyze Eq. (|3C|) for three particular cases. 

(i) If IV^ > r B * b is valid, then from Eq. (|l]) it fol- 
lows that TA'Ata =ln4, and Eq. ( |30| ) reduces to 

»i = \JC B * A \ 2 /T A » A . (32) 

(ii) For T a * a ~Tb*b the relation rs. B tp = 2 — \[2 
holds. Thus, Eq. (BG) reads as 



Wl ^|/c s , A | 2 2fV2-i) e -( 2 -^)/r 



(33) 



(hi) In the case where Tb*b ^^A'A is valid, one finds 
Tb'S^o = hi4, so that Eq. (30) reduces to 



\K-b*a\ 2 I^B'B ■ 



(34) 



For molecules, a rate regime is commonly considered, 
and the transfer rate is calculated by means of Fermi's 
golden rule in second-order perturbation theory with re- 
gard to the molecule-field interaction, 



w 



E 



PiWfi 



(35) 



where Wfi is proportional to the absolute square of 
the second-order interaction matrix element and the en- 
ergy conserving 5-function S(ujf — Wj), and the sum runs 
over the (quasi-)continuum of initial (i) and/or final 
(/) vibronic states. In particular, short-distance energy 
transfer, where the intermolecular coupling is essentially 
static, has been well known as Forster transfer Jl2|, O . 
Later on energy transfer over arbitrary distances has been 
considered (see, e.g., pi |l0| and references therein). 

For two two-level atoms, the single-transition probabil- 
ity per unit time Wif corresponds to P B . Starting from 
the effective Hamiltonian ([l3]) and evaluating Cs{t) from 
Eq. ( |l7| ) in first-order perturbation theory, we derive, on 
making the standard long-time assumption, 



Pb = 2tt \K,b-a\ 8(uja-£jb) 



(36) 



which is fully consistent with the result, obtained in 
second-order perturbation theory on the basis of the orig- 
inal (fundamental) Hamiltonian ([!]) The only new 
feature in the current treatment is that the medium in- 
duced atomic frequency shifts are taken into account. 

For establishing a rate regime, it is necessary that, ac- 
cording to Eq. (p5[), a continuum of initial and/or final 
states is involved in the transition. In the two-atom prob- 
lem at hand, the continua of states are obviously provided 
by the atomic level broadening due to the spontaneous 



decay. Assuming Lorentzian line shapes, we may perform 
the integrals in Eq. ( |35| ) to obtain 



4|£ 



B»A\ 



Fa*A + r_B*j 



■Pa(&a)- 



(37) 



Let us identify pa(u> a ) with Pa (to) [Eq. ( p7| ) together 
with Eq. (|3l])] and compare the result with the nonper- 
turbative rate w\ calculated from the exact temporal evo- 
lution of P B (t) [Eq. @ together with Eq. @]. Noting 
that Pa fa) ~ exp(— r^.^tp), we find for the three cases 
considered in Eqs. (p2) - (p4J) the following results. 



(i) r^- J 4>r B . B : 

Pi (to) 

(ii) T A , A =T B 'B- 
P A (t )^e-( 2 -^), 

(hi) T B *B^^A'A- 

P A (t ) 



w 



w 



V2- 1 ~ 0.41. 



1, ^ 

w 



(38) 



(39) 



(40) 



Agreement between w\ and w is observed in the first 
case, where Tb*b is sufficiently small. With increasing 
value of Tb* b an increasing discrepancy between w± and 
W is observed. The reason can be seen in the fact that 
there are two processes which simultaneously drive atom 
B: the energy transfer from A to B and the spontaneous 
decay. Hence, W\ is actually the rate of both processes 
combined, not the bare rate of energy transfer. This ex- 
plains why W\ is typically smaller than w and why the 
discrepancy between them becomes more substantial for 
enhanced spontaneous decay of atom B. To roughly com- 
pensate for the first effect, the ratio w\/w may be mul- 
tiplied by exp(TB*Bto). The values of the so corrected 
ratio are then ~ 1 for the cases (i) and (iii), and ~ 0.74 
for the case (ii). 



C. Strong atom-field coupling 

For the sake of transparency, we again consider identi- 
cal atoms that have equivalent positions and dipole ori- 
entations with respect to the material surroundings so 
that Eqs. J25| ) and ( p^ ) hold. Introducing the probabil- 
ity amplitudes 



C±(t) = 2-*[C A (t)±C B (t)] 



(41) 



of the superposition states |±) = 2~ 1 / 2 (\U A ) ± \Ub)), 
from Eqs. ([17]) we find that the equations for C+(i) and 
C_(t) decouple, 

C±(t) = dt' dLuK ± {t,t l ] Lu)e flSA ' Bi - t ~ t ' ) C ± {t'), 
Jo Jo 

(42) 
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where 



K± {t, t';u) = K A * A {t, t'; u) ± K A . B (t, t'; to) 



(43) 



Let us restrict our attention to the case when the ab- 
solute value of the two-atom term K A *s{t, t'; u>) is of the 
same order of magnitude as the absolute value of the 
single-atom term K A * A (t, t'; ui), so that there is a strong 
contrast in the magnitude of K + (t, t'; lo) and K-(t,t';uj). 
Assuming that the imaginary part of the Green tensor in 
the resonance region of strong RAFI has a Lorentzian 
shape, with uo m and Ato rn being the central frequency 
and the half width at half maximum respectively, we can 
perform the frequency integral in Eq. ( fl2| ) in a closed 
form in a similar way as in the single-atom case [[14], [l5| . 
In particular for exact resonance, i.e., to rn — Co a T &a*b, 
we derive 



C±{t) = 2-h- Au ™ t / 2 cos(n±t/2) 

(0± 3> Aw,„), where 

^± = x/2r±Aw m , r± = t a *a ± f a *b 



(44) 



(45) 



with T A * A and I 1 ^* b being defined according to Eq. pi] ) 
with uj m in place of Co a (—&b)- For the probability am- 
plitudes of the remaining states |=f), which are weakly 
coupled to the field, we obtain 



c T (t) = 2-h~ r ^ 2 

It then follows that 



(46) 



P 



MB) 



it) 



-rv* 



-Au) m t 



cos 2 (fi±t/2) 



+ (-) 2e- (A ^ +r ^ t/2 cos(fi±t/2) cos(26 A *Bt)] . (47) 

Note that in the equations given above the upper (lower) 
signs refer to the case where the state |+) (|— )) is strongly 
coupled to the medium-assisted field. 

For small damping, between three typical cases can be 
distinguished. 



(i) 4|fo. B |»n±, i<2/fi±: 

P A (t) = cos 2 (S A * B t), 
P B (t) = sin 2 (8 > B t). 

(ii) 4\6 A * B \~tt±, t4Zl/\2S A *B-n±/2\: 

P A (t) = \ [l + 3cos 2 (ft±i/2)] 
P B {t) = \ [l-cos 2 (ft±i/2)] . 



(hi) 4|5 A , B |<fi±, t<l/|2* A . B | 



cos 4 (ft±i/4) , 
sin 4 (f7±i/4) . 



(48) 
(49) 



(50) 
(51) 



(52) 
(53) 



From Eqs. (|48|) - (|53J) the following time-averaged proba- 
bilities P A , Pb, and P L = 1 — Pa — Pb are obtained, with 
the respective time integral being taken over one cycle, 
(i) Pa — Pb= \ and Pl=0. The excitation energy is 
periodically exchanged between the two atoms through 
virtual field excitations exclusively, (ii) P A = | , Pb = g , 
and Pl = § ■ The two exchange channels - one channel 
through virtual and the other one through real field ex- 
citations - compete with each other and destructively 
interfere, leading to a partial trapping of the excitation 
energy in atom A. (iii) P a = Pb= § , ■?£ = §• The in- 
teratom energy exchange is dominantly mediated by real 
field excitations. 



IV. POWER SPECTRUM 

Let us finally address the problem of the influence of 
the RDDI on the power spectrum of the light emitted 
during the decay-assisted energy exchange. To calculate 
the (physical) spectrum, we apply the formula (see, e.g., 

H) 



S(r,u s ,T) = / dt 2 / d*i 
Jo Jo 



-io)s(*2— ti) 



E(-'(r,t 2 )E(+)(r,t 1 )\l,(54) 



where ws is the setting frequency of the (ideal) spectral 
apparatus, T is the operating-time interval of the detec- 
tor, and 



E (+) (r) = / dcjE(r,w), 
Jo 



EM(r) = E< + !(r) 



A. Weak atom-field coupling 



(55) 



(56) 



For weak RAFI, we derive, on basing on Eqs. (JTs|) , 

ja, ©, and m, 



S(r,uj s ,T -> °o) = ^ 



Fa + F 



b 



Alo s + 5 A * B +iT+/2 
2 



' b 



where 



and {A' = A, B) 



ALu s -S A * B + iT-/2 
Aws = u s - Co A 



, (57) 



(58) 



neoc- 



/>'oo 

/ dwImG(r,r A sw)d^((u A »-w). 
Jo 



(59) 
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FIG. 1: The power spectrum 5[647r£oc 3 ro/(3fi(ii)] in the 
weak coupling regime for two atoms situated near a dielec- 
tric microsphere of single-resonance Drude-Lorentz-type [wt, 
transverse frequency; cup = 0.5u>t, plasma frequency; 7 = 
10~ 6 wr, absorption parameter; g( = 20At, sphere diameter 
(A t = 2ttc/ljt); u) m = 1.0504867 w T ; Aoj m ~5x 10" 1 lo t ; Ar A 
= ta — 0.5 d = Ars = 0.02 At, distance of the atoms from 
the sphere surface; cIa =ds, radially oriented real transition 
dipole moments; R — 20.04 At, interatomic distance; the de- 
tector is positioned at r — 20At on the radial line connecting 
the sphere center and atom A; uja = 1.05048 ojt (solid line), uja 
= 1.050485 w t (dashed line), and Q A = 1.05048621 w T (dotted 
line) . 



Note that replacing £(l>a' — (*A with it5(uja' —oS) would 
be too rough here. Equation ( p7| ) reveals that the RDDI 
may result in a doublet structure (two asymmetric lines 
at uja T Sa'B, with r± and \Fa ± Fb\ 2 being the widths 
and weights respectively). The line separation is seen 
to be twice the RDDI strength. A system of two two- 
level atoms, one of them initially excited, is obviously 
equivalent to a three-level system with two upper dressed 
states |±). The doublet structure can be understood as 
a result of the transitions of the dressed states to the 
ground state. 



As can be seen from Eq. (57), an experimental observa- 



tion of the doublet structure of the emitted light requires 
a delicate balancing act. The interatomic distance should 
not be too large to provide a reasonable level splitting, 
but it should not be too small to avoid |F^| = |Fs|, i.e., 
quenching of one of the two lines. It is worth noting that 
the presence of macroscopic bodies may facilitate the de- 
tection of the doublet, because it offers the possibility 
of realizing strong RDDI even for interatomic distances 
much larger than the wavelength. 

Another interesting feature is that, according to 
Eq. (p5|), either T + or T_ can be much smaller than 
^a*a(—^b*b)- That is, the resonant dipole-dipole in- 
teraction can give rise to an ultranarrow spectral line, 
albeit each time at the expense of the other line of the 
pair. If, e.g., a single atom is placed sufficiently near a 
microsphere, its spontaneous decay may be suppressed, 
with the emission line being accordingly narrowed. Com- 



pared to the emission line of a single atom, one line of 
the doublet observed for two atoms being present may 
be further narrowed by several orders of magnitude [ pT| . 

The RDDI-induced asymmetric splitting of the power 
spectrum is illustrated in Fig. |l] for two atoms near a 
microsphere. The solid line in Fig. |l| corresponds to an 
off-resonant atomic transition with not too strong RDDI 
and moderate contrast between T + and T-. As the 
atomic transition frequency uja approaches (from the off- 
resonance side) the frequency, where the absolute value 
of the RDDI strength reaches its maximum, the two spec- 
tral lines move apart from each other, with one of them 
gradually getting so much broaden that it eventually dis- 
appears (dashed and dotted lines). As uja shifts further 
towards the resonance line center, the RDDI weakens 
while the contrast between T + and T_ increases even 
more. As a result, the singlet structure persists with the 
peak moving back to the central position (not shown). 



B. Strong atom- field coupling 
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FIG. 2: The power spectrum S[64iceoc 3 r /(3huj A )} in the 
strong coupling regime for the three sets of parameters [To = 
10 _6 cjt, n+~128r ; (i) u>a = 1.04835747 u> T , i?~0.01A T , 
Sab - 2129 r (solid line); (ii) Cb A = 1.05045444 wt, R ~ 
0.027 At, Sab — — 32.2 To (dashed line); (iii) loa = u> m = 
1.0504867 wt, R= 20.04 A T , Sab ~ 0] (dotted line). Other 
parameters are the same as in Fig. |IJ 



For strong RAFI, Eq. ( |57| ) changes to, on basing on 
Eqs. @ and ©, 



S(r,uj s ,T -> 00) = - 



(W A ± W fl ) 
1 



Auj s ±6a-b+ n±/2 + iAu m /2 
Aojs ± 8a*b - n±/2 + iAoJm/2, 



8 



i(F A TF B ) 



where {A' = A, B) 



ALU S T5 A *B + iT T /2 



ImG(r,r A ',w m )dA» 



€qC 2 Q± 



(60) 



(61) 



Here the upper (lower) signs again refer to the case where 
the state |+) (|— )) is strongly coupled to the medium- 
assisted field. Eq. (60) reveals that due to the strong 
RAFI a doublet observed for weak atom-field coupling 
may become a triplet, with one of the lines of the dou- 
blet being split into two lines. These lines separated by 
Q± have equal widths (which are solely determined by 
the width of the medium-assisted field resonance) and 
equal weights. Note that their width and weight are dif- 
ferent from those of the third line, which is closely related 
to a line of the doublet observed for weak RAFI. From 
Eq. (^o|) it is also seen that, depending on the point of 
observation, this line or the strong-coupling-assisted dou- 
blet can be suppressed due to the interference effects. 

The effect of strong RAFI on the spectrum of the emit- 
ted light is illustrated in Fig. ^ for two atoms near a mi- 
crosphere. Three typical situations are shown, namely 
(i) stro ng, (ii ) moderate, and (iii) weak RDDI (cf. Sub- 
section [II C ) . Note that the triplet structure is most 



pronounced in the first case. 



V. CONCLUDING REMARKS 

We have developed a fully quantum mechanical the- 
ory to describe both RDDI and RAFI in the presence 
of arbitrarily configured dispersing and absorbing mate- 
rial bodies. In particular, the concept of RDDI works 
well for arbitrary atom- field coupling strengths. Though 
both the single-atom level shifts and the RDDI strengths 
result from the off-resonant atom-field interaction and 
are formally given by the similar-looking Eqs. (0) and 



dllf), respectively, Eq. (||) applies only to the level shifts 
caused by the scattering part of the Green tensor. The 
level shifts that are related to the vacuum field fluctu- 
ations have to be dealt with separately due to the di- 
vergence of the real part of the vacuum Green tensor 
at equal positions. In our calculations, they have been 
thought of as being already included in the bare atomic 
transition frequencies. By contrast, the real part of the 
Green tensor at different positions is finite and Eq. ( |l(i| ) 
can thus be used to calculate the total RDDI strengths 
caused by both vacuum and scattered-field fluctuations. 
In the same context, it might be incorrect to apply the 
rotating-wave approximation in the study of the RDDI. 
While the (scattering-assisted) level shifts b~\, A , which 
result from counterrotating terms, are small in general 
[cf. Eqs. (Q) and (||)], the corresponding RDDI strengths 
^A'A'(A^A') can b e significant, especially when the mu- 
tual distances of the dipoles are very small [cf. Eq. (|To|)]. 

To demonstrate the usefulness of the theory, we have 
applied it to the problem of resonant energy exchange 
between two atoms. In particular, we have shown that 
the RDDI can give rise to spectral line splittings of the 
light emitted during the decay-assisted energy exchange 
process. With a locally excited dipole as a probe, the 
effect may be employed, e.g., in imaging single atoms 
on a surface by scanning near-field optical microscopy 
techniques. The (numerical) calculations performed for 
atoms near a microsphere can of course be extended to 
other geometries of the material surroundings such as 
planarly stratified or cylindrical objects. It should also 
be mentioned that the theory is also readily extendable 
to the case of externally driven atoms. 
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